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Abstract 

We reformulate the second-order Schrodinger equation as a set of two cou- 
pled first order differential equations, a so-called "Shabat-Zakharov system", 
(sometimes called a "Zakharov-Shabat" system). There is considerable flexi- 
bility in this approach, and we emphasise the utility of introducing an "aux- 
iliary condition" or "gauge condition" that is used to cut down the degrees 
of freedom. Using this formalism, we derive the explicit (but formal) general 
solution to the Schrodinger equation. The general solution depends on three 
arbitrarily chosen functions, and a path-ordered exponential matrix. If one 
considers path ordering to be an "elementary" process, then this represents 
complete quadrature, albeit formal, of the second-order linear ODE. 

Keywords: Schrodinger equation, Shabat-Zakharov system. 
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1 Introduction 



The Schrodinger equation has now been part of mathematical physics for almost 85 
years [1]. Over the years it has been incorporated into many textbook discussions 
(for a necessarily selective subset see[2l[3lia[5l[6l[7l[8l[9l[T0l[ni[l2l[T3l[Tl[T^ 
it has been the subject of technical research monographs (for a necessarily selective 
subset see [HI [TBI [191 EDI EH [221 [23l El]), and considerable work has gone into 
exploring the mathematical foundations of the subject — see for example [251 [Ml 
[271 [28l [29] . Despite its long and venerable history, foundational questions related to 
the Schrodinger equation still periodically lead to new results [301 [SH [321 [321 [Ml [35] . 

In the current article we shall develop a very flexible formalism for reducing the 
second-order Schrodinger equation to a system of two first-order differential equa- 
tions, a Shabat-Zakharov system. This generalizes several of our earlier results, in 
particular those reported by one of the present authors in [30], and provides a formal- 
ism for formally solving the Schrodinger equation in terms of a 2 x 2 "path-ordered 
exponential" matrix. In addition to these formal developments, the technique is no- 
table for the fact that in appropriate circumstances it permits one to derive useful 
bounds on the behaviour of the exact wave- function. 



2 Basic idea 

Consider the one- dimensional time-independent Schrodinger equation [H [2, [3l [H [3, 
[6l[71[8l[9l[IOl[IIl[l2l[T3l[Tll[T5l[T6]: 



2m dx"^ 
Introduce the notation 



ft^ d^ 

tlj{x) + V{x) ^Ij{x) = E ^Ij{x). (2.1) 



2 2m[E - V{x)] 



HxY = ^ ^2 - (2.2) 



So we are really just trying to solve 

d^ 



dx^ 

or equivalent ly in the time domain 

d2 



^Pix) + k{xY ij{x) = 0, (2.3) 



ij{t)+uj{tyij{t)=0. (2.4) 



dt 
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Motivated by the JWKB approximation, 



, exp[z f k(x)] exp[— i f k(x)] 

the key idea is to re-write the second-order Schrodinger equation as a set of two 
coupled first-order hnear differential equations for the coefficients appearing in this 
linear combination. 

Systems of differential equations of this type are often referred to as Shabat- 
Zakharov systems [Sni EQl [35], or sometimes Zakharov-Shabat systems. A similar 
representation of the Schrodinger equation is briefly discussed by Peierls [21], and 
related representations are well-known, often being used without giving an explicit 
reference (see for example |36]). However an exhaustive search has not uncovered 
prior use of the particular representation presented here. (Apart, of course, from 
related precursor work by one of the current authors in [30], and our own more 
recent related work in ^Tl [331 [SH [SS]-) Nor, outside of our own work, has any 
attempt been made to use this Shabat-Zakharov representation to place rigorous 
bounds on the behaviour of one- dimensional scattering. 

We will start by introducing two arbitrary auxiliary functions (f{x) and A(x), 
which may at this stage be either real or complex, although we do demand that 
f'{x) 7^ 0, and then define: 

^(x) = aix) + bix) ^2.6) 



This representation effectively seeks to use quantities somewhat resembling the "phase 
integral" wavefunctions as a basis for the true wavefunction [2T]. We will ultimately 
want to interpret a{x) and b{x) as "position-dependent JWKB-like coefficients"; in a 
scattering problem they can be thought of as "position-dependent Bogoliubov coeffi- 
cients" . The representation given above is of course extremely highly redundant: one 
complex number ip{x) has been traded for two complex numbers a{x) and b{x), plus 
two essentially arbitrary auxiliary functions (p{x) and A(x). To reduce this freedom, 
or more precisely keep it firmly under control, we introduce an "auxiliary condition" 
(or "auxiliary constraint", or "gauge condition"): 

d /aexp(zA)^ ^ d / &expHA) \ ^ ^^^^ ^^.7) 



dx \ / dx \ -JTp 

Here x(2^) is yet a third arbitrary function of position. It is allowed to be complex, 
and may be zero. The original analysis, published in jBH] corresponds to the special 
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case A(x) = and x(^) = 0? so that it clear that the current analysis is a significant 
generalization. Subject to this "gauge condition" , it is easy to evaluate: 

^ = {a{x) exp{+i(f + iA) — b{x) exp(— — zA)} + x ^- (2-8) 
ax 

Repeated differentiation of this equation will soon lead to our desired result. 



3 Probability current 



To give us some insight into the physical meaning of the coefficients a{x) and b{x) it 
is useful to first calculate the probability current. As usual we take 



' 2m i \ dx dx ^ / m y dx^ 



(3.1) 



Here {h/m) is just a normalization (that is often set 1 for convenience). There is 
nothing really important in this normalization (unless we want to calculate experi- 
mental numbers), so we might as well set 

^{x,t) = lm(^r^y (3.2) 
Using our JWKB-based ansatz in terms of a{x) and b{x) we compute: 

^ = Im i^/^{a{x) exp{+iip + iA) — b{x) exp{—iip — iA)} + xip j-, (3.3) 
whence 

^ = Re < w — [a(x) exp(+i</? + iA) — 6(x) exp(— iy9 — -iA)] 

X [a{x)* exp{-iip* - iA*) + b{x)* exp{+iip* + iA*)] \ + lm{x} (3.4) 



This implies 

^ = Re 

+Im 




[|apRe{. 



,+i((p+A-<^*-A*) 



} - \b\^Re{e 



-i{(p+A-ip*-A*) 



}] 




Im {ab* e^(^+^+^*+^*)} + Im{x} 



(3.5) 



which we can finally recast as 



Re{/} n 
Im{(^'} 



a|2Re{e+2i'"(^+^)} - |6|2Re{e-2i'"(^+^)}] 




(3.6) 



Recall that at this stage v^(x), A(a;), and x{x) are completely arbitrary possibly 
complex functions subject only to the constraint 9?' 7^ 0. 

If we now temporarily demand that ^{x), and are real we see that 



an observation that strongly suggests that at least in those circumstances the quan- 
tities a{x) and h{x) might usefully be thought of as "position-dependent Bogoliubov 
coefficients" . 

4 Schrodinger equation as a first order system 

We shall now re-write the Schrodinger equation in terms of two coupled first-order 
differential equations for these position-dependent JWKB/Bogoliubov coefficients 
a{x) and h{x). To do this we evaluate the quantity d^-j/^/dx^ in two different ways, 
making repeated use of the gauge condition. From 



(3.7) 



ax \ vv'' 



d ( ■ ^' 




} + 



(4.1) 



dx^ 



we first see 




(4.2) 
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so that 



12 



+1^' { 2 



d /ae*^ 



But then 



dx^ 



dx \ ^/ip' 



2z(^' da . 



- /if' dx 



(4.3) 



+x 



(4.4) 



So finally 



iT^ dx 



(4.5) 



Now use the gauge condition to eliminate dajdx in favour of d6/dx. This permits 
us to write the quantity d^^/'/dx^ in either of the two equivalent forms 



d^^ 
dx^ 



[x' + x'-(^')']^-2^^ 



, dh e-^^-^^ 

dx 



1^' 



(4.6) 



and/or 



dV 
dx^ 



[x' + X - (V^')'] ^ + 2^V5' 



,da e+^^+^^ 
dx 



-2vVA'e''^+'^a - z [(^" + 2x<^'] h 



(4.7) 



We now insert these formulae into the Schrodinger equation written in the form 

d2^/; 



dx^ 



+ k{xy ^Ij = 0, 



(4.8) 



to deduce the first-order system: 



da 
dx 



+ 2^y [^'(^) + X' + X'- ivT - 2^' A'] a 



(4.9) 



dx 



-z[P(x) + x' + x'-(^')'-2^'A'] 4- 



(4.10) 



It is easy to verify that this first-order system is compatible with the "gauge condi- 
tion" (12.71) . and that by iterating the system twice (subject to this gauge condition) 
one recovers exactly the original Schrodinger equation. These equations hold for 
arbitrary v^(x), A(x), and x{x)i real or complex (subject only io ip' ^ Q to avoid 
divide by zero issues). 

This system can compactly be written in 2 x 2 matrix form, and once this has 
been done these equations exhibit a deep connection with the transfer matrix for- 
malism [37]. Let us define quantities pi(x) and p2{,x), not necessarily real, as 

(4.11) 



pi = i^" + 2xv'; 92 = k\x) + x' + X - ivT- 
We then re-write the Shabat-Zakharov system in 2 x 2 matrix form as 



d 

dx 



1 

V 



i[p2 - 2(p'A'] {pi + ip2] exp(-2z(p - 2zA)' 

{pi - ip2} exp(+2«(p + 2zA) -i[p2 - 2<p'A'] 



This has the formal solution 



a[x) 
b{x) 



E{x,Xo] 



a(xo) 
b{xo) 



12) 
(4.13) 
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in terms of a generalized position-dependent "transfer matrix" [37] 

E{x,Xo) = 



ex ^ 



da; 



(4.14) 



where the symbol V denotes "path ordering". 

Equivalently, if we were to be working in the time domain we would have 

E{t,to) = 



1[P2 - 2ip A] {pi + ip2} e-2*'^-2»A' 



dt 



(4.15) 

where T would now be the well-known "time ordering" operator (more usually en- 
countered in a quantum field theory setting) and we would now define 

Pi = (^ + 2x0; p2 = uj\t) + x^ -X-ivf, (4.16) 

with (p{t), A(t), and xit) being arbitrary functions of t rather than x, and 
k{x) — > u{t). 

Returning to position space, we can now write the (exact) wave function in inner 
product form 

il){x) = -^j== [exp(+i(y9 + iA); exp(— — iA)] 



a{x) 
h{x) 



(4.17) 



to yield a formal hut completely general solution for the Schrodinger equation 



iIj{x) 
Explicitly 

%Ij{x) = 



1 



= [exp{+i(^ + iA); exp{—i(p — iA)] E{x, xq) 



a{xo) 
b{xo) 



= [exp{+iip + iA); exp{—i(p + zA)] 



(4.18) 



(4.19) 



xPexp 



XO 



dx 



a(xoJ 
h{xo) 



i[p2-2^'A'] {pi + zp2} 
>i - ^P2} e+'^^+'^^ -z[p2-2v9'A'] 

This is the explicit general solution to the Schrodinger equation. It depends on the 
three arbitrarily chosen and possibly complex functions </?(x), A(x), and x{^)i ^-iid 
a path-ordered exponential matrix. If you consider path ordering to be an "elemen- 
tary" process, then this is a complete quadrature, albeit formal, of the Schrodinger 
equation, and implicitly, of the general second-order linear ODE. 



9 



5 Special cases 



We can now use the freedom in choosing (f{x), A(x), and xi^) to explore some special 
cases where the Shabat-Zakharov system simplifies. 



5.1 Case: A' = p2/{2ip') 

No one can prevent us from choosing 



A' 



P2 



that is 



which implies 



/\2 



A 



A 



, _ k\x) + + X' - iv') 

e{x) + x^ + x' -ivT 



dx. 



(5.1) 

(5.2) 
(5.3) 



Doing this greatly simplifies the Shabat-Zakharov system since now 
1 



da 
dx 

db 
dx 



That is 



^{ (b" + 2xif'] + I [k\x) + + x' - {^'f]) e-^^--^^^ ) b, 



(5.4) 



+ — { {[^" + 2x^'] - I [e{x) + + x' - (^')']) 



a. 



(5.5) 



d_ 

dx 



1 



{pi + ^^2} exp(— 2z(y9 — 2iA) 

{pi — ip-z} exp(+2z(y9 + 2zA) 





a 




b 



(5.6) 

In this situation one has eliminated the diagonal part of the metric — the cost 
unfortunately is that the off-diagonal components will now have a rapidly oscillating 
phase. (A somewhat related pair of equations can be found in reference |27j.) 
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5.2 Case: A = —Lp 

No one can prevent us from choosing 

A(x) = -V3(x), 

in which case 

1 



da 
dx 



_,[k\x) + x' + x'+{vT] &|- 
The comphcated phase structure has gone away, and we now have 



d 

dx 



1 



{pl - IP2] -l[p2 + 2(<^')'] 





a 




h 



(5.7) 



(5i 



(5.9) 



(5.10) 



The phases of the matrix entries are now slowly varying — the price to pay is that 
there is a full complement of nonzero matrix entries to deal with. 



5.3 Case: A = 

We include this case mainly for historical reasons, as it is an otherwise unpublished 
result that was the first significant generalization we obtained of the original result 
published in [30] ■ The Shabat-Zakharov system in this case simplifies to 



d 


a 


1 


ip2 


{pi + ip2}exp{-2iip) 




a 


dx 


b 


~ V 


_{pi -ip2}exp(- 


^2i(p) -i[p2 - 2ip] 




b 



(5.11) 



Since on still has the freedom to choose both (p and x this is definitely more general 
than [30], that article corresponding to the specialization x — > 0. 
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6 Application: 

Bounding the coefficients a(x) and b(x 

One of the particuarly interesting apphcations we have found for the Shabat-Zakharov 
system derived above is that it is possible to use it to place quite general and 
rigourous bounds on the coefficients evolution a{x) and b{x). From the general 
Shabat-Zakharov system 

(6.1) 

£ = +^|(b" + 2x^1-z[A:^(x)+x^ + x'-(^')^]) e+^^-+^^^a 

-z [k\x) + x' + x'- - 2^' A'] foj , (6.2) 

we see 

.da 1 



a 



, — 1^ [k\x) + x' + x'- {V'f - a* a 

+ {[^" + 2xv'] + ^ [k\x) + + X' - {V'Y] ) ^^^-^"^ oTh). 



(6.3) 



Therefore 

.da da* ^ k\x) + x' + x'-{v'f-'^v'^' , . 
a - — h a— — = Im < > a a 



+ Re 



dx dx I Lp' 

{[if" + 2xip'] + I [fc'(x) + + - {y'?]) e-^'^-^'^ g* h' 

(6.4) 
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But Re (A) < \A\, whence 



d|ap 
dx 



< Im 



+ 



Therefore 



d|a| 
da; 



< Im 



2^' 



(6.5) 



[^" + 2x^']+z[P(x) + xVx'-(v^'r] 



^2Im(</3+A) 



\h\. 



(6.6) 



While up to this stage A, and x have been allowed to be complex, we have found 
that for current purposes (establishing the bounds) it proves impractical to retain 
this level of generality, and to make any progress we must restrict attention to real 
A, and x- The inequality now reduces to 

[^" + 2xv']+i[k\x) + x^ + ^-{^'Y] 



d|a| 
dx 



< 



2^> 



and so 



d|a| ^ + 2xy? + +X' + X'- {y'fX 



\h\. 



(6.7) 



(6.8) 



dx ~ 2\^'\ 

We note that A has now completely disappeared from the inequality. Under the 
current assumptions it is easy to check that 

/ = Im(^/^V) = lap - 



so current conservation implies 



\a\^-\b\^ = l. 



(6.9) 



(6.10) 



(Ultimately, it is this equation that allows us to interpret a{x) and b{x) as "position- 
dependent Bogoliubov coefficients".) In view of this relation between a{x) and b{x) 



we have \b\ = y |ap — 1, so that we can deduce 

dial . \/ + ^XV? + [k'{x) + x' + X'- iv'n 



dx 



< 



2\ip'\ 



(6.11) 



13 



But this inequality can now be integrated. For convenience let us define 



^ = - r-, • 6.12 



Then 

d|a| 



dx 

But now 



<^?V|a|2-l. (6.13) 
dial 



v/|a|2-l dx 



dx< I d dx, (6.14) 



so that 



{cosh^^ |a| }^-^ < / dx. 

J Xi 



(6.15) 



Now apply suitable boundary conditions: as Xj — oo we can choose to set things 
up so that we have a pure transmitted wave, so |fe(— cxd)| = and |a(— (X))| = 1. On 
the other hand as xj — > +oo we must then choose to set things up so that a(x) and 
h{x) tend to a and /?, the Bogoliubov coefficients we are interested in calculating. 
Thus taking the double limit Xj — oo and xj ^ +oo we see: 

/+00 

dx. (6.16) 
'OO 

That is 

|a|< cosh ^9dx|. (6.17) 

This is the central result of this article — it can be modified and rearranged in 
a number of ways, and related inequalities can be derived under slightly different 
hypotheses, but all the applications we are interested in will reduce in one way or 
another to an application of this inequality or one of its close variants. 

For notational convenience, we often find it is useful to adopt the shorthand 

//• + 00 
= j , (6.18) 

since then 

\a\ < cosh < ® t9 dx > . (6.19) 
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From the normalization condition ( IG.lOp we immediately deduce 



T > sech^ 1^ dx| , (6.21) 
R < tanh^ i (p ^ dx\ . (6.22) 



\p\ < sinh 1^ i3 dxj . (6.20) 

When translated into equivalent statements about transmission and reflection prob- 
abilities, we flnd 



and 



where we reiterate 

ddx= i ^ — — dx. (6.23) 

I I 

The equivalent bound in the case x = was previously derived in [30], and via a 
rather different technique verifled in [31]. The current bound is deflnitely stronger 
than anything reported in [30l [31] , though somewhat surprisingly it can (after some 
transformations) be shown to be equivalent to the bound derived in [33] by using a 
radically different technique involving the Miller-Good transformation. Be that as 
it may, the underlying Shabat-Zakharov system is ultimately of deeper signiflcance 
and we continue to investigate the possibility of deriving improved bounds using the 
current and related techniques. 



7 Discussion 



There are several ways of recasting the Schrodinger equation in a form where it is 
more amenable to formal analysis. In this article we have rewritten the Schrodinger 
equation in terms of an equivalent system of first-order equations — a Shabat- 
Zakharov system — and then analytically studied this system. In particular we have 
used the system to derive a number of rigourous bounds on transmission probabil- 
ities (and reflection probabilities and Bogoliubov coefficients) for one-dimensional 
scattering problems, and compared them with earlier results in [301 l3Tl [33] . 

Even though the calculations we have presented are sometimes somewhat tedious, 
we feel that they are more than worth the effort — since there is a fundamental lesson 
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to be learnt from them. Technically, we demonstrated that the Schrodinger equa- 
tion can be written as a Shabat-Zakharov system, which can then be re-written in 
2x2 matrix form. We explicitly derived the general solution in terms of a position- 
dependent "transfer matrix" involving the symbol V which denotes "path ordering" . 
This explicit general solution to the Schrodinger equation depends on the three arbi- 
trarily chosen functions V5(x), A(x), and x(^) ^"^^ ^ path-ordered exponential matrix. 
If one considers path ordering to be an "elementary" process, then this is the holy 
grail of ODE theory (complete quadrature, albeit formal, of the second-order linear 
ODE). 
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